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Abstract. - Conformal field theory predicts finite-size scaling amplitudes of correlation lengths 
universally related to critical exponents on sphere-like, semi-finite systems S''"^ x R of arbitrary 
dimensionality d. Numerical studies have up to now been unable to validate this result due to 
the intricacies of lattice discretisation of such curved spaces. We present a cluster-update Monte 
Carlo study of the Ising model on a three-dimensional geometry using slightly irregular lattices 
that confirms the validity of a linear amplitude-exponent relation to high precision. 



Introduction. - The observation that the symmetry of systems at a critical point goes 
beyond scale invariance emerging from the divergence of correlation lengths, and additionally 
comprises rotational and translational as well as inversional invariance has greatly enlarged 
the scope of exact results that can be derived for such systems on field-theory grounds. While 
customary scaling theory and the real-space renormalisation ansatz |^] give rise to general 
relations between the scaling exponents and the concept of universality Q of exponents and 
certain amplitude ratios, conformal field theory (CFT) || allows for a classification of models 
according to their operator contents, thus providing the values of exponents and amplitude 
ratios, and, in special cases, even the amplitudes themselves. By interpreting the finite size of 
a system as an additional field in the scaling ansatz, these considerations apply not only to 
thermal scaling but as well to the scaling of observables of finite systems at a critical point in 
the limit of diverging system sizes, i.e. to finite-size scaling (FSS) Concerning universality 
in FSS, based on renormalisation-group considerations it has been argued that for semi-finite 
systems below their upper critical dimension the correlation lengths (measured in units of the 
lattice spacing) themselves are universal quantities Q. Exploiting conformal invariance this 
statement can be strengthened and the amplitude of this universal correlation-length scaling 
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can be calculated at least for a special geometry. To see this, consider the logarithmic map 

w ^ — Inz, z E C, (1) 
27r 

in two dimensions, which wraps the complex plane C around an infinite-length cylinder 
of circumference L. Exploiting the facts that the critical, connected two-point correlation 
function in the plane is fixed by the postulate of conformal invariance and that the given map is 
conformal, the critical two-point function of any conformally transforming (primary) operator 
on the cylinder 5^ x M follows immediately. In particular, correlations of the considered 
operator along the M-axis decline exponentially with a correlation length of 

where x denotes the scaling dimension of the operator. Thus, one arrives at a FSS relation 
including the amplitude A = l/27rx, implying universality of the correlation length itself. 

Since there is an essential difference between the conformal groups in two dimensions, 
where it happens to be infinite dimensional, and higher dimensions, where it is reduced to 
a finite-dimensional Lie group, very few CFT results exceeding those already known from 
renormalisation-group arguments survive a transition to higher dimensions. However, con- 
sidering the transformation (|l|) as a simple change of coordinates for a moment, rewriting 
it in polar coordinates allows for a generalisation to higher dimensions, now mapping R*^ to 
5'''"-'^ X R. Noting that such a map only acts on the radial part of the coordinates, leaving the 
angular part invariant, Cardy conjectured another linear amplitude-exponent relation for 
the correlation lengths in the geometry 5^"^ x R, namely: 

e„ - f . (3) 

where R denotes the radius of S"^^^ and x — x{d) is again the scaling dimension. Note that this 
includes the d — 2 result ^ with L = 2nR. This relation, if it can be confirmed, constitutes 
one of the few exact results for non-trivial three-dimensional systems. Since the notion of 
primarity of a scaling operator is a priori not well defined for d > 2, however, it is not entirely 
clear to which operators eq. (^) should apply. To find the analogues of primary operators 
for d > 2 one might start with the operator product expansion (OPE), cp. the explorative 
analysis of the operator content of the 0(n) models in three dimensions that can be found in 
ref. Note that this generalised mapping now connects different geometries, whereas the 
two-dimensional mapping can be understood as a meromorphism of the Riemann sphere onto 
itself. 

Considering alternative approaches, attempts of numerical investigation of the FSS of 
systems on such spherical geometries are severely hampered by the fact that there is only 
a finite number of triangulations, i.e. regular polyhedra, for each sphere S'^~^. A numerical 
transfer matrix study for the case of d = 3 using Platonic solids as a regular discretisation of the 
sphere was found to be inconclusive due to the restriction in system sizes 0. Here we follow a 
different line and use slightly irregular lattices as sphere discretisations for c? = 3, exploiting and 
explicitly probing for the universality of the considered observables. Considering universality, 
it might be favourable to study ratios of correlation lengths which are known to be quite 
insensitive to lattice distortions. The attempt to escape such complications by considering the 
flat geometry S'"'^ x 5^ x • • ■ x R instead surprisingly yields similar results, which are, however, 
up to now not connected to CFT calculations [[l2|. 



Lattice discretisation. - The most obvious model lattice with spherical topology is given 
by a cube covered by a rectangular mesh of lattice points, a lattice type that will be denoted 
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by (C) in the following. This inevitably introduces irregularities as compared to an idealized 
regular (large volume) triangulation of the sphere, which are given by a concentration of 
curvature around the cube corners and the defective coordination number of the corner points 
(three neighbours instead of four). In view of this problem several refinements have been 
proposed such as the replacement of the cube corners by triangular plaquettes, a construction 
that hides the coordination number defects away into the dual lattice, or the projection of the 
cube geometry (C) on the sphere by the application of appropriate weight factors to the 
links between lattice sites (S) This latter lattice can be shown to arrive at the "right" 

continuum field theory in the thermodynamic limit JT^ . On the other hand, even simpler 
approximations could be devised such as for example a "pillow" resulting from the glueing 
together of rectangles along their four sides. Since, at least for bulk quantities, differences 
between those different choices of model lattices appear to be small especially between 
(C) and its projected version (S), we here concentrate on the cubic discretisation (C) which 
is computationally much more convenient than (S). As far as universality is concerned, the 
ratios of asymptotic correlation lengths are found universal even with respect to marginal 
perturbations so that we do not expect to see artefacts of discretisation defects here. 
Amplitude universality, however, might be sensitive to such distortions. Other types of lattice 
discretisations will be considered in a forthcoming publication 

The cube geometry (C) is being realised by six L x L square lattices appropriately glued 
together along their edges. Since the model lattices including type (C) are only topologically 
equivalent to spheres one has to conceive a definition of effective radii to enable a proper FSS 
analysis in the sense of eq. (H). Assigning a unit volume to each lattice site, to each pair of 
bonds, or to each individual square of the model lattice, respectively, on counting these entities 
one arrives at areas of 

r 6L{L-2) + 8 "sites", 

1 6i(i-2) + 6 "bonds", "squares", 



from which the effective radii result via the relation R = \J AjA.'K. Since for a "quadrangula- 
tion" one obviously has 

4# (squares) = 2# (bonds), (5) 

the "bonds" and "squares" definitions coincide, but differ from the "sites" definition by a 
constant shift, leading to a slightly different approach to the leading FSS behaviour. 

Model and simulation. - We consider a classical, ferromagnetic, nearest-neighbour Ising 
model with Hamiltonian 



= - J ^ fJiCTj, (Ti = ±1. (6) 



The spins reside on lattices compound of the cubical sphere discretisation (C) with edge lengths 
L times a further linear lattice direction of length L^. We apply periodic boundary conditions 
in z-direction to eliminate surface effects in modelling the infinite R part. To control the effect 
of finite Lz still present we enforce the condition Lz/S,\\ 3> 1 in a systematic way, i.e. we scale 
1/2 proportionally to the radii R and thus (to leading order) proportionally to L. We use a ratio 
of Lz/£,\\ ~ 15 which empirically proves sufficient to let the effect of finite drop far below the 
threshold of statistical fiuctuations. Simulations were done at a quite precise recent estimate 
for the Ising model critical coupling in three dimensions /3c — 0.221 6544(3); the effect of 
critical coupling uncertainty was checked by a temperature reweighting technique and found 
negligible compared to statistical errors. For the FSS analysis we performed single-cluster 
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Monte Carlo simulations of nine systems with cube edge lengths ranging from L = 4 to 
L — 12, corresponding to effective radii of i? « 2 to i? « 8, and collected up to about 8 • 10® 
approximately independent measurements after an initial equilibration phase. Taking the 
scaled longitudinal sizes Lz into account, system sizes thus range up to about 3 • 10^ spins. 

Within the framework of CFT, the densities of magnetisation and energy turn out to 
be primary operators of spin models in two dimensions. Taking this as a hint for higher 
dimensions, we thus determine the longitudinal correlation lengths of these two observables 
to check whether Cardy's conjecture eq. (||) holds in the discretised case. We start measuring 
the corresponding connected, longitudinal correlation functions G'^'j'^(z), maximally enhancing 
the statistics by appropriate averaging and projection techniques, for details see [|l8j. We 
find that for the extraction of the correlation lengths from this information at a level of 
high precision the customary way of fitting the correlation function data to the functional 
form G(z)^^ll = a exp(— 2;/^||) + b for large distances z has severe drawbacks. Due to the 
statistical nature of the data one is not allowed to assume b = 0, i.e. the theoretical limit for 
infinite-length time series, a priori for finite-length measurements and is thus forced to use 
intrinsically unstable non-linear fits. Therefore, we resort to a diff'erence-ratio type estimator 
of the form 



■ G'^.ll(z)-(S"=-ll(z- A) 
G''^\\{z + A) -C^^Hz) 



(7) 



which eliminates the inconvenient constants a and b above. The parameter A > 1 controls the 
signal-noise ratio of the estimator and is being adapted for a minimum statistical error. To find 
final estimates for the correlation lengths ^|| (i?) of each system, we average the estimates ^|| (z) 
over a regime of distances z delimited by short distance deviations from the purely exponential 
behaviour and exploding fluctuations in the limit of very large distances. The averaging is done 
in a statistically optimized way entailing an expensive analysis of variances and covariances as 
well as autocorrelations, based on the "jackknife" resampling technique JTsf . 

Results. - Collecting the final estimates for the correlation lengths £,\\,a/e of the densities 
of magnetisation and energy we arrive at the finite-size results shown in fig. |^(a), which scale 
linearly to leading order in the effective radii R defined above. Differences between the two 
radii definitions can hardly be distinguished at this scale. Plotting the effective amplitudes 
£,\\/R, however, reveals the presence of corrections to the leading scaling behaviour that are 
clearly resolvable at the given level of accuracy, cp. fig. 1(b). We thus fit our data to a FSS 
ansatz including a leading correction term, 

^\\{R)=AR + BR", (8) 

treating a as an additional fit parameter that thus constitutes an effective correction exponent 
including parts of the higher-order corrections as well. This higher-order effect is kept small 
by successively removing points from the small R end while monitoring the quality-of-fit 
parameter Q. The range of points used is indicated by the range of the fit lines in fig. ^ 
Combining the fits to the data for the two different radii definitions we arrive at final estimates 
for the leading correlation length scaling amplitudes and their ratio of 

= 1.996(20), 

A, = 0.710(38), (9) 
A^/A, = 2.81(15). 

Using weighted averages of recent estimates for the critical exponents of the three-dimensional 
Ising model fll], namely v = 0.63005(18) and 7 = 1.23717(28), the conjectured amplitudes for 
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Fig. 1. - (a) FSS plot for the spin-spin correlation length „.(/?). (b) Scaling of the amplitudes 
^W^rr/R- The continuous lines show fits described in the text and the horizontal line indicates the 
conjectured amplitude A^J"^' = 1.9298(13). 



comparison are 

A^'ij = 1/xa = 2/{d - -ijv) = 1.9298(13), 

A^'^j = l/a;, = l/(d - \lv) = 0.70780(23), (10) 
^conj/^conj ^ ^^j^^ ^ 2.7264(13), 

where d = 3 is understood, cp. eq. (^). Comparing the amplitude ratios, which are known to 
be more robust with respect to lattice inhomogeneities ||l^, |l^ , we find very good agreement 
between the simulation and Cardy's conjecture. The same holds true for the amplitude of 
the energy-energy correlation length whose estimate is naturally less precise than that of 
the spin-spin correlation length. For the latter we still obtain reasonable agreement, but 
with a deviation of the numerical estimate towards larger values which is on the edge of 
statistical significance. To check whether this result indicates a systematic deviation due to 
the irregularity of the used lattice discretisation (C), simulations for the other discretisations 
mentioned above in the introduction should be performed | p6[ . 

Conclusions. - Approximating the sphere S'^ by a cube covered by a rectangular mesh of 
lattice points, we find that Cardy's conjecture of a linear amplitude-exponent relation in the 
FSS of correlation lengths on the sphere-like geometry 3"^^^ x R holds true for the case of the 
three-dimensional Ising model and the densities of magnetisation and energy. In terms of the 
analogy to CFT in two dimensions the considered operators are thus analogues of "primary" 
operators in three dimensions. The ratio of leading scaling amplitudes does, at the given level 
of accuracy, not show any effect of the irregularity of lattice discretisation; for the amplitude of 
the spin-spin correlation length further simulations using different discretisation schemes are 
necessary to judge the extent of universality that is being obeyed. Our results indicate that 
the geometry 5*^ x K considered here is closer to the original cylinder 5^ x R (as far as the 
prevalence of universal features of correlation lengths scaling is concerned) than the toroidal 
geometry x R which is not conformally flat, but nevertheless exhibits a scaling law quite 
similar to that found in two dimensions 
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